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ABSTRACT 


The single-blow transient testing facility of the Naval Postgraduate 
School was reconstructed. With this apparatus the flow friction and 
heat transfer characteristics of four grades of steel wool, considered 
to be representative of fibrous metals, were determined. The flow 
friction characteristics were determined on the basis of the isothermal 
pressure drop across the medium, while the heat transfer results were 
determined by means of the single-blow transient testing technique. 

The heat transfer results are presented as the Colburn j-factor. 

The flow friction results are presented in two forms: as the Fanning 


friction factor, f; and as the correlating friction factor, C All 


£g 
results are presented as a function of Reynolds number based on the 


product of the experimentally determined permeability and tortuosity 


parameter of each sample. 
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NOMENCLATURE 


English Letter Symbols 


a 


Linear coefficient of Reynold's 
equation, reciprocal of the 
permeability (1/K) 

Matrix total heat transfer area 
Matrix minimum free flow area 


Mateix total frontal- area 


Matrix solid cross-sectional area 
available for longitudinal conduction 


Quadratic coefficient of Reynold's 
equation, interpreted as a tortuosity 


parameter 


Cubic coefficient of the empirical 
equation for the pressure gradient 


Correlating friction factor 
Filtered 00 heat at constant pressure 
Matrix material specific heat 
Matrix thermal capacity (Wc) 
Rod diameter in stacked screen model 
Fiber hydraulic diameter 
Flow passage hydraulic diameter 
Friction power per unit area 
Mass flux (m/A ) or (m/A. ) 

e fr 
Gravitational acceleration 


Proportionality factor in Newton's 
Second Law 


Surface heat transfer coefficient for 
convection; heat transfer power per 
unit area per degree temperature dif- 
ference 


Typical Units 


F7 


dimensionless 


dimensionless 
BTU/(1lbm-°F) 
BTU/(1bm-°F) 
BUF 

ft 

ft 

fit 

hp/ft* 
1bm/(hr-£t”) 
UN 


32.2 حا‎ p s 


BTU/(hr-ft“-°F) 





Fluid thermal conductivity 
Matrix permeability 
Contraction coefficient 
Expansion coefficient 
Matrix thermal conductivity 
Length parameter in Noo 
Total matrix flow length 

Mass flow rate 

Pressure 

Matrix porosity na cp 

Heat transfer rate 

Specific gas constant (53.35 for air) 
Flow passage hydraulic radius (D,,/4) 
Temperature 

Absolute temperature 

Normalized temperature 

Flow velocity 

Bulk velocity (6/0) 

Specific volume > 

Matrix volume 

Average fiber width 

Fluid mass in matrix 

Matrix mass 

Uncertainty interval for quantity y 


Distance along flow passage from the 
matrix inlet 


Pitch of rods in stacked screens model 


Reduced length No, 0/13) 


BTU/(hr-ft-°F) 
gc? 
dimensionless 
dimensionless 
BTU/(hr-ft-?F) 
ft 

fe 

1bm/hr 
lbf/ft* 
dimensionless 
BTU/hr 
ft-lbf/(lbm-9R) 
Et 

SE 

SR 
dimensionless 
ft/sec 

ft/sec 

Fe an 

p 

in 

1bm 

1bm 


units of y 


"E 


dimensionless 


dimensionless 





Greek Letter Symbols 


o = € CG Se 


Specific surface (compactness)(A/A, /L) gel 
Ratio of orifice diameter to pipe dimensionless 
diameter 
Difference or change (of time, tem- 
perature, etc.) 
Time sec, hr 
Fluid viscosity lbm/(hr-ft) 

| 1bf-hr/ft? 
Free time dimensionless 
Density EE 
Kinematic viscosity TT 


Subscripts 
atm Local atmosphere 
ave Average 
f Fluid (gas, air) 
1 Initial; inlet 
m Matrix; mean 
O At orifice 
S Solid (matrix material) 
STD Standard (temperature and pressure) 
x Local conditions 
1 Inlet conditions (upstream of matrix and heaters) 
2 Inlet conditions at matrix entrance 
3 Exit conditions at matrix outlet 


Dimensionless Groupings 


E 


Correlating friction factor; ratio of dissipative 
effects to inertial effects 





Pr 
Re 


St 


Se 


Fanning friction factor; ratio of wall shear stress 
to fluid dynamic head 


28 


Colburn j factor (No N_ ). This factor plotted versus 


Lo 
Reynolds number defines the surface heat transfer char- 
acteristics 


Longitudinal heat conduction parameter for solid material 
0 9 
Time parameter (hA9/W_c_) 


Prandtl number; ratio of viscous dissipation to thermal 
dissipation (i.e., kinematic viscosity/thermal diffusivity) 


Reynolds number; ratio of inertial effects to viscous 


یں سیب 


Stanton number; ratio of heat transfer power per unit 
surface area to fluid thermal capacity per unit flow 
area (h/Gc ) 


Number of transfer units; ratio of total heat transfer 
power to total fluid thermal capacity Dire c 
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T. INTRODUCTION: 


Regenerators add on the order of ren percent to the thermal effi- 
ciency of the basic gas turbine power cycle. Due to the current 
emphasis on the use of gas turbines, as well as the reawakening interest 
in the Stirling and Ericsson Bes, where a regenerator is an integral 
part of the cycle, much attention has been paid to the design of these 
heat exchangers. Frequently the applications of such cycles are weight 
or volume-limited, and therefore regenerators must possess a high heat 
transfer surface density. 

However, as noted by Kays and London in their monograph on compact 
heat exchangers po) increasing the heat transfer area necessarily 
increases the pressure drop (hence the pump work required) across the 
heat exchanger, particularly when the working fluid is a gas. Thus, 
the consideration of the heat transfer characteristics of a surface for 
compact heat exchanger applications cannot be divorced from a consid- 
eration of the friction characteristics of the surface. 

Between 1965 and 1968 several compact heat exchanger surfaces were 
extensively studied at the Naval Postgraduate School. These included 
perforated nickel plate matrices and offset rectangular fin matrices 
of the AiResearch Manufacturing Company studied'by Piersall [23]; six 
different plate-and-fin cores by the Solar Division of International 
Harvester studied by Ball [2], Traister [29], and Trost [30]; and the 
ceramic cores by Corning Glass, studied by Bruno [5]. All of these 


سے 


E Numbers in brackets refer to items in the bibliography. 
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surfaces were fundamentally porous to flow in one direction only and 
were intended for use in gas turbine plants. 

Fibrous metals (such as steel wool) possess the high surface density 
requisite for a compact heat exchanger. Moreover, the apparent small 
hydraulic diameter of the effective flow channels of such materials 
should lead to large values of the convective heat transfer coefficient 
h, since h always varies with a uus power of the channel hydraulic 
Bieneter [13]. Also the tortuous nature of the three-dimensional flow 
through fibrous media hinders the development of the momentum and thermal 
boundary layers, which should further increase the value of h. Although 
the high surface area to volume ratio of the individual fibers excludes 
fibrous metals from high temperature environments (e.g., the open cycle 
gas turbine regenerator), such materials do have application in low- 
temperature closed power cycles and in refrigeration cycles such as the 
reversed Brayton, reversed Stirling, Gifford-McMahon, or Pulse Tube 
C los. 

bis tesis was directed at four objectives: 


1. Reconstruct the transient test apparatus (dismantled since 
1968), including new instrumentation. 


2. Test a previously studied compact heat exchanger surface 
to insure reliability of the reconstructed apparatus. 


3. Determine the heat transfer and flow friction character- 
istics of various grades of fibrous metals, specifically 
steel wool. (The "grade" of a fibrous metal refers to 
the size of the fiber, and not necessarily to the porosity 
of the fibrous metal matrix. Porosity was not an inde- 
pendent test variable). 


4. Ascertain the appropriate length parameter to be used in 
the definition of a Reynolds number upon which a unified 
correlation of the heat transfer and flow friction char- 
acteristics may be based. 
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A. REYNOLDS NUMBERS FOR POROUS MEDIA FLOWS 
Flow friction and heat transfer data are usually correlated with 
the flow Reynolds number for the surface geometry under consideration: 
جو ان‎ pale 
where 1 is a characteristic length parameter for the surface geometry. 
The results thus correlated are independent of any particular geomet- 
rically similar surface. For flow in fibrous media the dilemma is the 
selection of the characteristic length parameter. If the fibers were 
packed so tightly that the flow could be unambiguously considered 
internal in nature, the obvious choice would be some sort of effective 
hydraulic diameter of a representative flow channel. This was the 
procedure in the previously tested cores (e.g., the Solar series, 
Cercor, etc.), which were porous to flow in one dimension only and 
possessed a well-defined channel hydraulic diameter. At the other 
extreme, 1f only one fiber of the medium were under consideration 
(i.e., purely external flow), the characteristic length parameter 
would be the hydraulic diameter of the fiber. However, neither such 
simple case exists. 

Three approaches were made to the selection of the characteristic 
length parameter. In the first, fibrous metals were assumed to be 
porous to flow in only one dimension and the equations applied to the 
previously tested cores were used. The second approach is based on 
the in-line stacking of crossed-rod matrices. Both of these idealized 
models produce a physically conceivable characteristic length. The 
third approach was based on the fundamentals of flow in porous media, 


but produces no such conceivable dimension. 
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1. Uni-dimensional Porosity Model 
. Consider a heat exchanger with n uniform channels of length L. 
LE A. is the free flow area of the entire cross-section, and A is the 


heat transfer area of the entire exchanger, then: 


Free flow area/channel A /n 
Nominal volume/channel - (A /n)L 


Interior surface/channel = A/n 


But the interior surface of a channel equals the cross-sectional perim- 
eter times the length of the channel. Therefore 

Perimeter/channel = A/(nL) 
The hydraulic diameter of the channel is defined as four times the cross- 


sectional area of the channel divided by the wetted perimeter, or: 


- 4A) , ¿Al 
i (A / [nC ]) A 


For such a model which is porous to flow in only one dimension, the 


(1) 


porosity (i.e., volume of voids over nominal volume) is: 
(P A.) Ae (2) 


where Ber equals the total frontal area of the medium. Since the 


nominal volume of the exchanger is the product of A " and L, the 


B 


specific surface (i.e., heat transfer area per unit of volume) is, 


with substitution from (1) and (2): 


= > oe id as 3 
Ó AL AD, A (3) 


r n 
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As explained in Appendix B, the porosity and heat transfer area (hence 
also the specific surface) may be determined independently of these 
relations for fibrous metals. Thus the ‘equivalent channel hydraulic 
diameter for this model may be computed from (3). This is the charac- 
teristic length used in the Reynolds number under this idealization. 
2. In-line Crossed Rod Matrices Model 

This model consists of square-mesh cross-rod screens which are 
stacked so that there 1s no separation ا‎ layers and so that the 
mesh in one layer coincides with the meshes in adjacent layers (see 
inset). The rods are assumed to be of diameter d and separated by 
uniform distance x, d. 

On the basis of a repeatable unit cell (similar to the procedure 


in the preceding section) the following relations are easily derived [13]: 


p1- Tx, © 
A. /A, 7 C مب‎ 
B= 1/4) de 
Des dx sd (6) 





The pitch, x,, is calculated from the known porosity, ^ with the use 


L. 


of (4). Again this equivalent channel hydraulic diameter is used as the 
length parameter for this model. Note that this model differs from the 
previous in that it is porous to flow in three dimensions, rather than 
one. For the cases where the fiber is not of circular cross-section, 
the rod diameter d is replaced by the fiber hydraulic diameter in the 


equation for the specific surface (5). 
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3. Flow in Porous Media 

Analysis of the Prenat tals of flow in porous media provides 
a length parameter which is more justifiable, in the sense that no 
geometrical approximations are made. 

In 1856 peer. after a study of the flow of water through the 
sand filters of the water purification plant at Dijon, France, made the 
empirical observation that the bulk (i.e., nominal) flow rate of fluid 
through the filters is proportional to the difference between the fluid 
inlet and outlet — and inversely proportional to the thickness of 
the filter. Since this observation was valid only macroscopically, 
and therefore over finite distances, it was later formulated on a dif- 
ferential basis as what is now known as Darcy's Law of fluid flow in 


porous media: 


D= -ROP A9) | (7) 


Here U is the bulk velocity of the fluid flow (i.e., based on the total 
Frontal area of the medium and not the internal free flow area), p is 
the fluid pressure in the medium, and g is the gravitational accelera- 
tion, The constant K is known as the permeability of the medium. 
Darcy's Law is the fundamental relation describing the flow 
of fluids in porous media. The significance of this law is that it 
relates the internal pressure field to the bulk velocity U, a quantity 
which is measured external to the medium. Together with the continuity 


equation for fluid flow in a porous medium: 


^ 3 (0? : v (aut (8) 


(uP = total porosity of the medium). 
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Darcy's Law completely specifies the flow field. Under steady state 
conditions these two equations then reduce to Laplace's equation, and 
a solution is possible. It must be emphasized that Darcy's Law is an 
empirical relation; it has no analytical basis. Theoretically, if a 
satisfactory model can be formulated for the ene medium, and suf- 
ficient boundary conditions stated, the same result could be achieved 
by a.solution (closed form or numerical) of the Navier-Stokes equations 
(which would involve the internal velocities rather than the external, 
bulk velocity). E practice this analytical approach has rarely been 
successful, largely because of difficulties with the boundary conditions. 
It is convenient at this point to define some terms of impor- 
tance in the study of Darcy's Law: 


a) Porosity,/2, is the volume of voids divided by the 
total (i.e., nominal) volume. 


b) Permeability, K is, a measure of the ease with which 
a fluid may transverse a medium under the influence 
of a given pressure gradient. 11 is important to 
note that the permeability has units of length squared 
and is independent of the fluid.  Permeability is 
defined by Darcy's Law (7). 
c) Tortuosity is the relative average length of the 
flow path of a fluid particle from one side of the 
porous medium to the other. 
Clearly the permeability is the important parameter of the 
“medium for predicting the pressure drop across it. Unfortunately, it 
is the porosity of the medium, and not the permeability, which is 
generally known. For a given porous medium, the permeability is a 
unique, but unknown function of the porosity [22, 27]. Several theories 
have been advanced which purport to theoretically relate porosity and 


permeability. The oldest is that due to Kozeny (1927) but his own 


experimental data shows discrepancies from the theory of -68.7 to +85.7 
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per cent [22]. The most recent theory is based on statistical mechanics, 
but regretably attaches the randomness to the fluid flow, and not the 
medium, where it properly belongs [27]. Thus it may be seen that the. 
most accurate determination of the permeability is by experiment. 

If Darcy's Law (7) is written for one-dimensional flow in the 
absence of significant changes in a potential, it becomes 


(after rearrangement): 


0 لام‎ . (9) 
dx K 


Obviously the permeability, K, can be determined by measuring the iso- 
thermal pressure drop across a thin slice of medium at different bulk 
velocities, U. (Constant temperature is required so ےج‎ remains 
constant). 

Mother re has not been so kind, however, for Darcy's Law is 
strictly valid only for low flow rates. At higher velocities the 
magnitude of the pressure gradient increases with velocity more rapidly 
than the linear rate indicated by Darcy's Law. Based on experimental 


data gathered in the non-Darcy regime, Reynolds (1900) suggested the 


relation: 
dP. au ل2‎ (10) 
dx 7n Jc 


Forchheimer (1901) added a cubic term to better fit the data. However, 
the coefficient of the cubic term was found to be of such small magnitude 
that the term is generally neglected. Missbach (1937) proposed that the 
pressure gradient in non-Darcy flow be proportional to U to some exponent 
between one and two. As this relation could not even be بج‎ 


tified theoretically, it has generally been ignored in favor of Reynold's 


0 








equation. The most refined formulation is that of Polubarinova-Kochina, 


which includes momentum changes due to compressibility effects: 


2 
E um. BM 2adl (11) 
dx y dc Je dO 
This additional term is the acceleration across the medium due to the 
decrease in the density of the fluid. Since both a T frame of 
reference and a frame moving at the I velocity are inertial frames, 


this time-dependent acceleration may be rewritten as a convective accel- 


eration [8]: 


| FIWM bol” dU 
Fr ai U + Te T £ s (11b) 


Although (10) is of experimental origin, it was initially jus- 
tified on quasi-theoretical grounds by analogy to the pressure drop in 
pipe flow. In laminar pipe flow the pressure gradient is proportional 
to the velocity (the Hagen-Poiseuille equation), while in turbulent pipe 
. flow the gradient is proportional to the square of the velocity. Since 
transition between laminar and turbulent flow in fibrous media would not 
be as sharply defined as in pipe flow, the reasoning went, the pressure 
loss must be the sum of the two analagous losses [22]. 

It will be recalled that the non-linear (viz., inertial) terms 
in the Navier-Stokes equations do not become significant until the - 
Reynolds number exceeds some critical value [28]. This critical value 
marks the end of the "slow flow" regime. Separation of the boundary 


layer in external flow over bluff bodies and similar effects in curved 
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channel flow occur in this new regime. Por both internal and external 
flows, there is a dramatic increase in the pressure drop long before 
the flow becomes turbulent itself. It is now thought that non-Darcy 
flow is due to the emerging influence of these inertial terms [9, 26]. 
Hubbert [12] proposed such a completely laminar theory wherein the 
increasingly important inertial forces (which increase with velocity 
more rapidly than do the viscous forces) distort the streamlines of 
flow. This leads to increased velocity gradients and therefore, in- 
creased negative pressure gradients (by Bernoulli's Law). These 
increased velocities also lead to increased viscous forces. The effects 
are lumped as inertial resistance, hence a greater pressure drop than 
that predicted by Darcy's Law. 


Dimensional analysis of Reynolds equation 
| ur ~~ 
-FẸ - ابره‎ + (10) 

shows that a must have dimensions of one over length squared MUI B while 
D mas divenovons otvenewover length (1/1) Inwfaet; forssmellebullsves 
locity, U, the quadratic (i.e., inertial) term in Reynolds equation 
becomes negligible, and (10) reduces to the one-dimensional Darcy's Law 
(9). Therefore the constant a must be one over the permeability (1/X). 
Thus (10) becomes: 

d boU 

Jul, bl 

E - K oe O. (12) 
2۳۰۴۰۰۰7۳۰7۳761۲ ط‎ of the inertial term is an inertial resistance coet- 


ficient thought to be ude to the successive contractions and expansions 


of the flow in porous media. Green and Duwez [9] assert that b is 
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actually a measure of the tortuosity of the medium. However, there is 
no purely theoretical reason for assuming that b is independent of the 
rend [3]. 


If (10) is arranged in the following form: 


l [d j- 4 1 

" 7 U | ye Tyg tbla: (13) 
it is seen that a plot of pu 7 ae.) versus (Ea. 2 ( should be linear 
with slope equal to the tortuosity parameter, b, and intercept equal to 
one over the permeability, 1/K. Beavers and Sparrow [3] have examined 
the 9۷ت‎ of water through a fibrous metal product known as FOAMETAL, 
which is structurally similar but not identical to steel wool. Their 
investigation covered a velocity range of 20:1 for 5 samples and showed 
no deviation in excess of 2.5% from the linear plot predicted by (13). 
This is considered to a dT indication of the validity of Reynolds 
equation and suggests a simple means of determining the permeability and 
tortuosity parameter for a porous medium. 

However, the use of this technique requires the approximation of 
ae by pfi, (in other words, a uniform gradient between pressure 
taps is assumed). This approximation may be eliminated by the integra- 
tion of Reynolds equation (10) or the Polubarinova-Kochina equation (llb) 
under the assumption of isothermal flow of a perfect gas. Since (11b) is 
the more general of the two, it is used for the illustration of the inte- 


gration procedure. 


ge sau AE bel + pe کا‎ (11b) 
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The viscosity is the absolute viscosity in units of force-time/length 


squared. Multiply (11b) through by the density 


pe =p + De 


(14) 


But = PART for a perfect gas, and لام‎ - ۵ : the mass flux. Substitution 


of these identities in (14) yields: 


Bn =apG + EC ts (Er) رجا‎ (15) 


The mass flux, G, is constant across the medium since it is based on the 


frontal area, m 


MT ge ayu E + + e. (E) (RT) (16) 


(Tuc. OT m/A. ). 02, 





- E =P= - تایبه.‎ + eG + E P RSE (17) 
O G + BG + S p(-1/p*) $2 (18) 
- cP = aul + BG = = (4 de) (19) 
EG +7 8x 2 ($e) 0 


Equation (20) can be integrated from entry to exit of the medium (a length 


D): 


23 








2 = 
2 2 3 
-PojaopGl + با هط‎ . E In (p) ۱ (21) 


l 





(Ga 2)_ G | [B )=a Gier b "x. ( 
ZRT q EL FF = 

This equation may also be arranged so that the slope is again equal to 

the tortuosity parameter b while the intercept is one over the permea- 


bility (a = 1/K): 





A |p-e _G In (Bi) A ba )23( 
ns Ze Je Pe jh 
The choice between (13) or (23) to determine the permeability 
depends on the conditions of the experiment. The former is restricted 
to isothermal flow with a constant pressure gradient between pressure 
measurement stations. The latter is restricted to isothermal flow of 
a perfect gas, but places no further restriction on the nature of the 
pressure gradient. 
Once the permeability has been determined for a medium, a 
Reynolds number can be defined. Since Nro is the ratio of inertial 
effects to viscous effects, and since all the viscous effects are 


embodied in the linear term in (10) or (11), while all the inertial 


effects are in the quadratic term: 


N _ L > J E - U Kb) (24) 








pus was previously identified as the absolute viscosity in force 
pros pg is the absolute viscosity in mass units. To avoid ambiguity, 


substitute the kinematic viscosity > TP . Therefore, 


U 

NL > pw (Kb) =: (25) 
It is seen that the characteristic length parameter based on this ap- 
proach is the product of the permeability and the tortuosity parameter 
of the medium. 

There ۳ obvious disadvantages to this whole procedure in 
defining a Reynolds number. First, a length parameter involving the 
product of the permeability and some vague tortuosity parameter cannot 
be conceived as'a physically identifiable length. Second, the medium 
must undergo a series of experiments before this pseudo-length parameter, 
Kb, can be determined. Third, even the most subtle difference between 
the internal geometries of two cores can make a unified correlation 
between them based on a Reynolds number defined by (25) impossible. For 
example, Beavers and Sparrow found that the presence of loose fiber ends 
within a medium made impossible any correlation with a medium that did 
not have loose fiber ends. 

Nevertheless, the proper length parameter is the one which cor- 
relates friction and heat transfer data into single curves, and this is 


the ultimate test.. 
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E. FRICTION FACTOR 

The friction factor is the dimensionless variable which correlates 
the pressure loss behavior of similar surfaces. It is loosely defined 
as the ratio of dissipative forces to inertial forces. Any force (and 
not exclusively viscous forces) which causes a pressure drop across a 
medium is a dissipative force. 

Again there is more than one approach. The first is the solution 
for the familiar Fanning friction factor, f, based on a channel analysis. 
The second is a friction factor which is an outgrowth of the definition 
of the Reynolds number in terms of the permeability, Cee 
۲ ۳۳۱۳۳ Friction Factor, f: 


Consider the following idealized situation [13]: 
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The pressure drop at entrance consists of two parts. First 
there is a pressure drop due to the decrease in the flow area, which 


from continuity results in a flow acceleration. 


لم - ll]‏ کے ج P‏ 


Note that the identity (A /A c) equals the porosity (2) is strictly 
true only for the model which is porous to flow in only one direction. 
For the other two approaches this identity slightly overestimates E 
since for these approaches the medium is porous to flow in three dimen- 
sions. 

Second, there is a pressure loss caused by boundary layer 
separation, due to the irreversibilities of an abrupt contraction, and 


manifested by a vena contracta. 


G 
Op = K. £— E = Ko 229e (27) 


where K is the empirical contraction coefficient. The total pressure 


loss at entrance is therefore: 





(ap) = E eS £ + Ke) (28) 


Similarly there is a pressure change at exit consisting of two 


parts. First, there is a pressure rise due to the increase in area: 


-ap = -A 1-p*) = 





23.08) (29) 
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Second there is a pressure loss due to the irreversibilities of a 


sudden expansion 
2 


UE 
Op“ e P e zag. (30) 


where K. is the empirical expansion coefficient. The total pressure 





change at exit is the sum of (29) and (30): 


(OP), ex Wa +e) (31) 


1 


1 





E aed «n? 
EE بو‎ 
(Kas l+ 
Pn Je 
The remainder of the pressure drop across the heat exchanger is 


due to the overall fluid acceleration and the core friction. The pres- 


sure loss due to overall fluid acceleration is given by: 


du‏ ہوا 
E Tees dx‏ 


C 0090۳0٢ CONST., 
جو زد‎ Sx Cho ) (33) 


(32) 


dx 


Mp)... = G* Cle, (Pa) (34) 


The pressure loss due to core friction defines the friction factor: 


u f(A 
an Ra. Ar 


(35) 


G P(e (36) 
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After replacement of density by specific volume, the total 


pressure drop across the heat exchanger is: 


^ 


(K.+1-e:) TUA 1 ln - (4- KIS (3‏ ی 
)2 
Zoe ` | e 8 e‏ (۵۳) 
For the isothermal flow of a perfect gas, the pressure is inversely‏ 
proportional to the specific volume and (37) may be solved for f in‏ 


terms of pressures: 


—— E + KN (BIEL. ak eco 


For heat exchangers where there are successive contractions and 
expansions (such as fibrous materials) the entrance and exit coefficients 
are effectively lumped into the core friction factor [13]. Therefore 


in this study: 


I 29 pl an Brel Ë E 1+ YE (39) 


For small pressure differentials, the first term in (39) is by 
far the greatest contributor to the friction factor. A further approxi- 
mation then is: 


PE = Pm = P, = Pa 


f = [290 Pm (E + P. «ehe 


(40) 


Equation (40) defines the Fanning friction factor for fluid flow in porous 


fibrous media. 
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2. Correlating Friction'Factor, Ce: 
Since the friction factor is the ratio of dissipative to inertial 
forces, the common practice in porous media studies is to divide the 


pressure gradient as given by Reynolds equation by one-half of the 
inertial term and define the quotient as a correlating friction factor 
I 9j: م١‎ 
C dx 
= cc mapa 
f لا‎ 4] 
2 bp (41) 
: y Ge 
Recall that the Reynolds number in terms of the permeability and tortu- 


osity parameter is: 


— = C (Kb) (25) 
Substitution of Reynolds equation 


23 | (12) 


into (H1) yields 


Zee | al LU 
ارط سای‎ "c | (42) 


)43( یہی 


The kinematic viscosity is re ~ Thus 


3 
meinl 


e سے‎ < == Va | (44) 














2 
It should be noted that لام‎ 3 Y and that therefore 


equation (41) for C, might be written .(with ep replaced by AP/L ) 


C, = 329: pom rs 3 ws (45) 


The term in braces is the first term in brackets in the expression for 
the Fanning friction factor, f, equation (40). Therefore the correlating 


friction factor is related to the Fanning factor by: 
A 
3 els AP (1+ *) ies 


CS HEAT TRANSFER 

The parameters which characterize the ability of a surface to transfer 
heat include the heat transfer power per unit heat transfer area (h), and 
two non-dimensional forms of h: the number of transfer units 0 and 
the; Colburn j-factor. The heat transfer power per unit heat transfer 


area, h, is actually the experimental constant in Newton's Law of Cooling: 


` = hA (AT) (47) 


and is more commonly termed the "convective heat transfer coefficient." 


The number of transfer units, N, , is the ratio of heat transferred by 


tu 


the surface per unit of time to the thermal capacity rate of the fluid: 


= 48 
N,, = hA mep) (48) 
The Colburn j-factor is a heat transfer parameter based on Reynolds 
analogy and empirical observation which is found to behave in a nearly 
identical manner as the friction factor. for most fully-developed internal 


flows [6]: 
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2 2/3 
J P NN Dir. , 3 NE N, (49) 

Due to instrumentation difficulties inherent with complex geometries, 
the usual steady state method of determining these characteristic param- 
eters are ea, een, for compact heat exchanger surfaces, 
and it is necessary to resort is transient یمم‎ In this "first- 
cut" ES of ee from fibrous metals, the simplest such 
eee tue. the "single blow transient test" or "maximum slope" technique 
Was used. This method has been amply described in references [2, 5, 11, 
155 16, 73 19, 23, 24, 25, 29, 30, and 31]. 

The problem of analytically determining the fluid and solid temper- 
atures as a function of position and time for a porous solid with constant 
initial temperature and step change in the fluid entering temperature was 
studied independently by Anzelius, Nusselt, Hausen, and Schumann between 
1926 and 1929. The Anzelius solution is overly simplified, while 
Schumann's solution is valid only for liquids. 

The approach of Hausen, however, appears to be valid for the problem 
at hand. If the fluid is flowing under the influence of a pressure gra- 
dient in one direction only, the problem may be lumped in the other two 
coordinate directions, and a one-dimensional energy balance may be made 


on a differential element of the solid. 


h G ts) malte + Hear) 
سس‎ 


ncp te 
(38x Zn, =ks As r sts dx) 
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Assumptions made are: 


a) 
b) 
c) 


d) 


e) 


Properties of the fluid are temperature-independent. 
Fluid flow is steady. 
Porous solid is homogeneous. 


Thermal conductivities of both fluid and solid are infinite 
perpendicular to the fluid flow direction. 


Thermal conductivity of fluid is zero in the flow direction. 


The boundary conditions are: 


a) 


b) 


c) 


The matrix is initially at uniform temperature. 


At time equals zero there is a step change in the entering 
fluid temperature. 


Matrix bowndaries are adiabatic. 


The energy terms (on a per unit time basis) are then: 


a) 
b) 
c) 


d) 


Ac, (S) 
Heat absorbed by solid =P ss 2 X 
Heat transferred to solid by convection - hit, E Bel 


at 
Heat transferred from the fluid by convection - me (S£ Jax 





. l | at 
Heat transferred from solid by conduction = - k A. (S ax 


Therefore energy balances on the solid and fluid are: 


RASSE Ske dy = kA + بل(‎ Es) dx 


(50) 


eps L .)dx = O (51)‏ ۳ که ره 


The following non-dimensional variables are now introduced: 


a) 


b) 


Reduced length, £ = ^ (| = Nu p 
Generalized time, T = € É 1 ہے‎ 





_hA Q ۳ HA Co X 
TW 00 ۳ 
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itudinal ducti = k, As 
c) Longitudinal conduction parameter, À= FE T 
Since the thermal capacity of the fluid (ee) within the solid matrix 
at any instant of time is much less than the thermal capacity of the 


solid, the second term in the definition of is negligible compared to 


the first. mM e ure. 
Ts ie, ^ م‎ (52) 


In terms of Du Z, and T equations (50) and (51) become: 


B | = ts 





(53) 


= - 2 er (54) 


For the case of zero ae conduction (i.e., A= 0), the above 


equations become simply: 


dts I 
TT Escrits uh 


n 


9t 
az 


This system of differential equations, subject to the initial and 


UAE | (56) 


boundary conditions previously specified may be solved analytically 
using Bessel functions. Theoretical fluid temperature - time curves 
at constant Z (in particular at the matrix exit, where Z^ Nor e Zn) 


may be plotted. To determine the N u of a heat exchanger surface, the 


experimental temperature-time history at the matrix exit is matched to 
one of the family of theoretical curves. This procedure is not only 
tedious, but also depends on the assumption of zero longitudinal con- 


duction. 
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In 1950 Locke [19] eliminated the necessity of matching experi- 


mental and theoretical curves. A new time variable is defined: 


J - T/N,, (57) 


He differentiated the analytical solution to Hausen's equations (retain- 


ing the assumption A= 0) and evaluated the derivative at the matrix exit 


(Z=N_): 
tu 
-N (1+ 
E ba ") 

s Li DU D - J, 2iN 3 e (58) 

ل 
where te = fluid exit temperature E 0, te - fluid exit temper-‏ 

2 5 1: 


N 


Ture for < 0, and m - steady state temperature. Locke showed that 
as long as the change in the fluid entry temperature is a step change 
(and longitudinal conduction is zero), the u of a heat exchanger sur- 
face is uniquely related to the maximum slope of the temperature-time 
curve (i.e., the maximum value of (58)) at the matrix exit, hence the 
name "maximum slope technique." 

Although the assumption of zero longitudinal conduction is essen- 
tially valid at large flow rates (i.e., high Reynolds numbers), it is 
quite invalid at low Reynolds numbers for most surfaces. If these 
longitudinal conduction effects are included in the analysis, the 
maximum slope is a function of À as well as ال‎ With finite differ- 
ence techniques, Howard [11] has included these effects and has tabu- 
lated (Table II) and plotted (Figs. 1,2) maximum slope as a function of 
N, o. However, for fibrous metals such as steel wool, it is felt 
that any possible longitudinal conduction path would be too lengthy 
for there 2 be any pronounced effect. Thus in this study Howard's 


results are used only for the special case of A = Or 
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Kohlmayr [15, 17] showed that the maximum slope technique may be 


4 02:9; and is singular 


inaccurate for 3.02 N, 7 2.0, unstable for N. 


at M = 2.0 (Fig. 3). He proposed an se curve matching technique 
based on the first moments of the areas under the experimental and theo- 
retical temperature-time curves to circumvent this difficulty [16].  More- 
over, this "centroid technique" is not restricted to step ce in the 
fluid entry temperature. As will "2 in the Experimental Tech- 
niques — in this study the step change is approximated by securing 
the electrical power to the heaters in the flow stream. The actual 
response of the air is exponential rather than stepped and therefore 
Kohlmayr's centroid technique could be used to correct this deviation 
from Locke's idealizations. However, Bruno [5] has studied the various 
transient testing techniques and reports that there is no measurable 
improvement in the results by using the centroid techniques for و‎ 
Due to the expected large surface area and restricted flow rate through 
fibrous metals, it was felt that all Non values would exceed 3.5, there- 
by eliminating the need to correct for deviation from step. 

Thus the heat transfer parameters New? j, h for fibrous metals were 
determined in this study by Locke's maximum slope technique, subject to 
the assumptions of exact step change in the fluid temperature, and zero 


longitudinal conduction. 
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11۳۰ DESCRIPTION -OF “SURFACES 


Steel wool was selected as a representative fibrous metal. Steel 
wool is low carbon bessemer wire with high tensile strength. It may 
be a continuous piece up to 100,000 feet long. It usually has three 
m. but may have five [4]. The federal government has set speci- 
fications on the fiber size and type of steel for seven grades between 


0000 and 3. The restrictions on fiber size are listed in Table I [7]. 


TABLE ۰ E requirements for individual fibers 


Not more than 5 Mean width Not more than 5 No fiber 

Class percent under of fiber percent over to exceed 
| Inch  — Inch Inch Inch 
0000 0.0006—0.0010 0.003 0.006 
000 .0008— .0013 .004 .008 
00 .0012— .0022 .005 .010 
0 0. 0003 .0016— .0024 .007 .012 
1 0.0005 .0018— .0026 .010 .015 
2 0.001 | .0040— .0060 . 014 025 
3 .002 .0070— .0095 .018 | .040 


a sana 


Grades 0000 through. 0 must be low carbon steel. Grades 1 through 3 may 
be low carbon steel or AIST 430 stainless steel. 

Grades 00, C, 1, and 3 were selected for testing. In the case 
of Grades 1 and 3 the type of steel was determined by a simple test. 
The plain carbon steel will take the copper out of a CuSO, solution, 
while A1S1 430 stainless steel will not. All grades tested were plain 
carbon steel. Sample fibers of the four grades were examined under a 
microscope with a movable platform. With the cross-hairs in the eye- 
piece and the platform movement control dial (graduated in hundredths 
of a millimeter) the fibers were checked against the federal specifi- 
cations. Cross-sections of each fiber wens examined to determine if 
the samples were three or five sided. All samples indeed were three 


sided. In particular, the samples of Grades 00 and 1 (which were 
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manufactured by the Brillo Corporation) appeared to have a very nearly 
equilateral cross-section. | 

Such an equilateral cross-section was assumed for all grades tested, 
with side length equal to the average of: the range of mean widths of 
the fiber quoted in the federal s esas ss (Table 1). The volume 
of the material in each sample is equal to the weight of the sample 
divided by the density of the steel. The assumed cross-section and 
volume of steel specify the heat transfer ۸ A. (See Appendix B). 

The steel Wool was tested at essentially the same porosity as that 
with which it is packaged. The porosity may be crudely estimated from 
article 5.1.1 of [7] which states that one pound of steel wool shall be 
packaged in a sleeve.approximately four inches in diameter and fourteen 
to sixteen inches long. Based on a density for plain carbon steel of 
490 lb Ft Sl, the porosity should be between 0.979 and 0.983. The 
porosity may be more precisely determined from the weight of the sample 
tested and the volume of the test cell in which it was packed. The free 
flow area is approximately determined from the porosity. 

Complete geometrical descriptions of each sample, based on the 
assumptions stated, are found in Figs. 13, 15, 17, € 19. 

A core known as Solar No. 4 tested by Bruno [5] was also tested to 


insure the reliability of the reconstructed apparatus. 
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IV. EXPERIMENTAL TECHNIQUES 


The single blow transient test apparatus (Figs. 4,5) was designed 
to conform to the assumptions made for the analytical solution of the 
heat transfer problem (equations (50) and (51)). These idealizations 
are: 

l. The velocity and temperature profiles of the fluid 

flow at any cross-section in the matrix are both 
steady and uniform. 

2. The thermal conductivity of the matrix is finite 

(or zero) in the direction of the fluid flow and 


infinite in tie direction normal to the flow. 


3. The thermal capacity of the matrix is much greater 
than that of the fluid contained within it. 


i. The thermal prepertaessef the fluid and matrix 
are both constant and uniform. 


5. Phe convective heat transfer coefficient is some 
suitable average and remains constant. 


6. A step change in the inlet fluid temperature occurs 
at real time equals zero. 


A specially designed nozzle, flow-straightening screens, and an even 
distribution of heater wires across the channel insure that the first 
assumption (i.e., uniform velocity and temperature profiles) is satisfied. 
Piersall [23] made velocity and temperature traverses at the matrix inlet 
section over the entire flow range and concluded that the velocity 
profiles were indeed uniform and that the temperature profile varied 
by no more than + 0.5°F. Wheeler [31] has demonstrated the importance 
of the uniformity of these profiles. 

For fibrous metals it is assumed that there is no effective solid 


area, E available for longitudinal conduction. This has the same 
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effect as assuming zero matrix conductivity in the flow direction and 
therefore e: the second assumption. 

As long as the fluid is restricted to ھ0‎ the assumption of 
negligible fluid thermal capacity compared to matrix thermal capacity 
is SOIT It must be noted that the testing procedures and results 
of this study are not applicable to liquids. 

By limiting the temperature rise to 209F the assumption of constant 
properties is met. The worst offender is the viscosity, which may vary 
D as + 1.5 percent from its value at the mean temperature. The 
limited temperature rise also insures that the convective coefficient, 
h, is effectively constant. 

As has been discussed in Section II, the assumption of a step change 
in the fluid inlet temperature cannot be precisely satisfied. In the 
current testing method the step change is simulated by turning off the 
heaters. This actually leads to an exponential decay in the air temper- 
ature. Techniques are avallable to compensate for this deviation from 
step change, but Bruno [5] has indicated that these corrections are not 
required for نا‎ 3.5. Since it was expected that the N. range 
experienced would exceed this value, no effort was made to so compensate. 

The data required for the computation of the various heat transfer 
and flow friction parameters is as follows: 


- atmospheric pressure (inches Hg) 


atm 
۲ orifice static pressure (inches H,O) 
OP - pressure drop across orifice (inches H,O) 
OP, - pressure drop across matrix (inches H,O) 
P - static pressure at matrix inlet (inches H,0) 
is fluid temperature at orifice (millivolts) 
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d. - orifice diameter (inches) 


= 


P - patio of orifice diameter to pipe diameter 
CS - chart speed (millimeters/second) 
t,-t, - downstream temperature response (mm) 


Pressures were measured with water manometers and checked with 
inclined gages where the flow rates and orifice plates permitted. An 
exception was the atmospheric pressure, which was measured with a mer- 
cury barometer. 

Only the matrix downstream temperature response, (t,-t,)> was 


measured, although the upstream temperature (t ) was monitored to 


Zu 
check for the 20%F temperature rise from the heaters. A more complete 
description of the equipment and instrumentation is furnished in 
Appendix A. A sample output is shown in Fig. 6. 

Because the turbocompressor which induces the draft through the 
apparatus, would not function for much more than seven minutes with- 
out overheating and tripping the circuit breaker, and because of the 
importance of isothermal flow to the permeability studies, the friction 
factor data runs were made separately from the heat transfer data runs. 
This was in contrast to previous studies where the friction factor data 
was collected immediately after the heat transfer data. 

For a heat transfer data run the pressure drop across the ASME 
standard orifice plate corresponding to the desired mass flow rate 
was predetermined. The number of heaters required to achieve a 20°F 
temperature rise at that flow rate was also predetermined. When the 
desired pressure.drop had been obtained, the recorder pens were zeroed, 
and the heaters were energized. The (t,-t,) channel indicated imme- 


diately whether the 20°F rise has occurred. When the core had reached 
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a steady state temperature (as indicated by the (t,-t,) trace) all 
pressure measurements and the orifice temperature were recorded. The 
heaters were then de-energized and the transient temperature response 
downstream of the matrix (t,7t,) was recorded. After the completion 
of the run the سے‎ slope was visually determined with a straight 
edge from the recorder output and entered in the data sheet. The data 
sheet conformed to the input format for a digital computer program 
written by Trost [30] and modified for fibrous ات‎ which reduced 
the data. (See Appendix B). 

A friction factor data run has made in the same manner as a heat 


transfer data run except that no heaters and therefore no recorder 


were required. 


42 











V. PRESENTATION OF RESULTS 


For each of the four grades of steel wool ec the heat transfer 
and flow friction characteristics have been computed and appear in 
Tables III to X. Although Fanning friction factors were computed 
for the heat transfer data runs (Tables VII to x), the flow state at 
the time of measurement cannot be guaranteed to have been isothermal 
(particularly at low flow rates). For this reason the friction factor 
results in Tables 111 to VI are preferred, since these runs were made 
without the heaters and are the basis for the permeability analysis. 
The Colburn j factor is plotted as a function of Npe for — sample. 
The values of پت‎ and h as a function of Nee appear in Tables VII to X. 

Figure 7 is a plot of the friction data in the form of (23). Figure 
8 is the same plot in the form of (13). From these plots the sample 
permeabilities and tortuosity parameters are determined. This informa- 
tion is prerequisite for the selection of a length parameter for the 
Reynolds number and therefore appears first in the results. Figures 9, 
10, and 11 plot the Fanning friction factor, f, as a function of Nee 
based on'the three different approaches attempted. The correlating 


friction factor, C „, is plotted against Nro based on Kb: in Figure 12. 


f 


Figures 13, 15, 17, and 19 contain the geometric and physical char- 
acteristics of the four samples. Figures 14, 16, 18, and 20 graphically 
present the heat transfer and flow friction characteristics (i.e., j and 
f versus N ۳ based on Kb) for each grade of steel wool. Figure 21 


R 


compares the j values of the various cores as a function of N 2 based 


R 
on Kb. Figure 22 compares the heat transfer power per unit of area of 


each core as a function of the friction power per unit of area. 
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Vi. DISCUSSION OF RESULTS 


A. RELIABILITY OF APPARATUS 

The Solar No. 4 core previously tested by Bruno was tested in the 
E 20 to insure its reliability. An unsatisfactory 
amount of electrical noise was present in the thermocouple signals, 
apparently due to a combination of the unshielded 0.001 inch thermo- 
couples (Appendix A), ultra-sensitive Brush amplifier, Variac power 
supply, and Physics Department linear accelerator, located immediately 
below the apparatus. With this system, Nog values were found to deviate 
by as much as 23 percent and generally 12 percent from Bruno's plot. 
The installation of a more sensitive reference junction (see Appendix A) 
may be part of the explanation. Friction factor results deviated by 
as much as 12 percent and were consistently 7 to 8 percent higher than 
Bruno's plot. 

Siieemeic Uncertainties were 1 1 TO Iepereent In سال‎ values Tmi? 


percent in N o values, and + 3 percent in f values for Bruno [2, 5] and 


R 
somewhat higher in the reconstructed system due to the previously 
mentioned difficulties (see section on experimental uncertainties), it 
was felt that these results constituted satisfactory agreement. In 
other words, the problem was not felt to be with the thermocouples, 
and the research proceeded. 

Subsequently a less sensitive amplifier-recorder system was installed, 
the Variac was replaced by a motor-generator DC source, and testing was 


performed when the linear accelerator was not in operation. Also a more 


elaborate manometer system was installed (Appendix A). 
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B. PERMEABILITY ANALYSIS 

The permeabilities and tortuosity parameters for the four samples 
were determined by means of a first order polynomial least squares curve 
through the data plotted in the forms suggested by both (13) and (23). 
The former equation assumes تا‎ 20 7 IT. flow with a con- 
stant pressure ne This technique was used by Beavers and Sparrow 
in their study of water flow through bore fibrous media [3]. The 
latter equation assumes isothermal flow of a perfect ase but is not 
restricted to incompressible flow or constant pressure gradients. 6 
permeabilities computed by the second method are on the order of one 
percent higher than those computed by the first method for the four 
samples examined. Therefore, compressibility effects are considered 
to be negligible in the computation of permeabilities and tortuosity 
parameters. This is not unexpected, since the maximum change in the 
pressure across the cores tested did nct exceed 50 inches of water. 

In Fig. 7 the results of the second type of analysis (i.e., com- 
pressibility effects included) are plotted. With the possible exception 
of grade l, a straight line plot is a reasonable curve through the data 
points. Outside of the data points corresponding to very small flow 
rates (where the relative error in the pressure measurements is large) 
and those corresponding to very large flow rates (where the absolute 
error is large, due to fluctuations in the flow), the maximum deviation 
of the data points from the linear least squares plot is ll percent. 
This ^ not nearly as good as the 2.5 percent maximum deviation reported 
by Beavers m TM [3]. However, their study was for the flow 
of an incompressible fluid (water) over a velocity span of 20:1, with 


14 pressure taps, cathetometer, calming reservoirs, and direct flow 
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measurement. In this study the fluid was compressible, the velocity 
span was P ensi; 200:1, with indirect flow measurement (orifice plates), 
2 taps, water manometers read by eye, and no calming reservoirs. With 
all this taken into consideration, Fig. 7 might be considered to be in 
adequate agreement with the empirical Reynolds equation (10). 

Deme ly. the data points show a marked concavity rather than 
random deviation from the linear least squares plot. This would indi- 
cate the presence of higher order terms than those in (10). Recall 
that Forchheimer proposed the addition of a cubic term to Reynolds 
equation, but that the coefficient was usually so small that the term 
was negligible except at very large velocities. In light of this and 
the extensive velocity range encountered in this study, a second order 
polynomial Meese Squares curve was passed through the data points plot- 
ted according to the manner suggested by equation (13). (Hereafter 
compressibility effects are ignored.) The results appear in Fig. 8 
along with the linear least squares curve. In general the data points 
deviate from the, second order least squares curve by less than half 
the amount of deviation from the linear curve. This is consider to 
be strong justification for the replacement of Reynolds equation (10) by 
Forchheimer's equation: 

U” 2 
-Jp yU bel +U (59) 
dx K Je 
The second order polynomial least squares curves plotted in Fig. 8 are 
representations of (59) recast in the following form: 


s 


TE Da Eyn + 
fa IY 32 (PU) E x (60) 
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c. (61) 

I MN CH b a? emm c 
ER K 3 a (62) 
(63) 


4152) ndis Gime) 


where compressibility effects have been neglected and 


E ce" fee p (64) 
= cp | (ge? (65) 


* 


It is not clear what physical significance this additional term has, 
but it evidently involves a coupling of inertial and viscous effects 
and is probably dependent on the fluid. Dimensional analysis shows 
that c* has units of (force-time 1 ine ) and that c is dimensionless. 


The particular results from Fig. 8 are: 


Grade Permeability Tortuosity 

Steel K 2 Parameter, b 

Wool (EE) AFE) c 

با- = 

00 ea x 1C 416.8 -9.57 x 10, 
0 G Ol 243.1 -3.02 x 10", 
1 SO KU 201.4 -3.98 x 10, 
3 6-517 x 0 61 -0.512x 10 


From these results the characteristic length parameters (Kb) are 
0 6000199 ft ١0.0000167 ft., 0.0000302 ft., and 0.00009"6 PE. for 


Grades 00, 0, 1, and 3 respectively. 
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The presence of this cubic term raises questions regarding the 
derivation of the Reynolds number based on the product of the permea- 
bility and a tortuosity parameter (25) where one term was said to embody 
all the viscous effects while the other embodied all the inertial effects. 
It is now evident that there is a third term which may involve both 
effects. Nevertheless, the whole procedure was an a DES one, 
and if the resulting Reynolds number correlates the data there is no 


point in questioning the validity of the derivation. 


C. CORRELATING LENGTH PARAMETER AND FRICTION FACTOR 

The Fanning friction factor, f, was chosen as the test for the 
proper length parameter to be used in Reynolds number. Once the hydro- 
dynamic correlation is: obtained, it was thought that the heat transfer 
correlation would involve the same length parameter. In Fig. 9 a cor- 


` 


relation of f and No. is attempted based on a flow passage hydraulic 
diameter from the model porous to flow in one direction. In Fig. 10 
the same correlation is attempted based on the crossed-rod matrices 
model. Figure 11 is the same correlation based on the product Kb. 

It is evident that the second model is better than the first in 
that the curves are somewhat closer together. However the plot of the 
curves f versus Neg based on Kb represents the tightest n e jug 
particular note that the curves for Grades 00 and 0 practically coincide. 
These Es grades had practically the same characteristic length (1.99 x 
107^ IE aNd 1.67 X o ft. respectively). The same is true of the 
curves for Grades 1 and 3 which also have comparable characteristic 
Vengtns (3.02 x ia ft. and 3.46 x 10^ Ft. -— ہ‫‎ This is 
considered to be a clear indication of the superior correlation to be 
had by using a characteristic length equal to the product of the medium 


permeability, K, and the medium tortuosity parameter, b. 
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The fact that a perfect correlation was not obtained in Fig. 11 
may be partly explained by the inverse proportionality between f and 
the heat ECT area, A. This area has the largest uncertainty of 
any variable in the study (see section on experimental uncertainties). 

The correlating friction en Cs is plotted against Noo 111107 


12. As was demonstrated in the section on theory, 1f Reynolds equation 


(10) is valid then C, is automatically determined by the Reynolds number 


f 
based on Kb: 


C, ۲ IN "de (44) 


The data points show adequate agreement with (44) for Reynolds numbers 
less than one. However, for Reynolds numbers greater than unity, equa- 
tion (44) is seen to over-estimate Ce as defined by equation (51). This 
is no surprise, since there appears to be a cubic term in the empirical 
equation for the pressure gradient. This cubic term is insignificant 


except at large velocities, which explains the deviation of the data 


points from the theoretical curve in Fig. 12 at high Reynolds numbers. 


D. HEAT TRANSFER AND FLOW FRICTION CHARACTERISTICS 

The heat transfer (4) and flow friction (f) characteristics of 
Gvades 00, 0, 1, and 3 are presented in Figs. 15, 16, 18, and 20 respec- 
tively and in Tables III to X. Table VII shows that for Grade 00 when 
j<0.00675, the corresponding value of Hes is less than 3.0, where the 
maximum slope method is highly unreliable. Similarly for Grade 0 when 
j<0.0118 and Grade 3 when j< 0.0579. All heat transfer results which 
fall in ENS noes must be disregarded, due to the high uncertainty 


in the-N values. 
tu 
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Also the, (t,-t,) Tec MII bur the last five sets of Grade 00 


data had an excessive amount of noise, and the determination of the 


maximum slope was at best a crude process. Inasmuch as the "noise" 

died off as the flow rate o cde TI is thought that it is actually 
the random fluctuations in the — similar to the fluctuations 
observed when the test section was removed and only the quiescent ambient 
air temperature was recorded. This ne perhaps indicative of some flow 
disorder downstream of the Grade 00 sample, but there is insufficient 
information to draw any firm conclusion. Thus all but the data points 
taken at the five largest flow rates for Grade 00 must be regarded as 
qualitative only. 

The Colburn n is the result of an empirical correction to 
the Reynolds analogy for non-unit Prandtl numbers [6]. The Reynolds 
analogy assumes that for fully developed flow (particularly turbulent 
internal flows), the hydrodynamic and thermal boundary layers are 
Similar. Therefore the heat flux at a surface is proportional to the 
shear stress at that same surface for a given fluid [8].. Thus the 
Stanton number is proportional to the friction factor and the curves 
of f and j as functions of Reynolds number are parallel. 

In previous studies of compact nes exchanger surfaces, the flow 
passages were well defined channels, and the parallel behavior of the 
f and ee predicted by Reynolds analogy (with Colburn's correction) 
was clearly observed [29]. This was not found to be the case for any 
of the four samples tested (Figs. 14, 16, 18, € 20). Whereas the 
friction factor was found to decrease with increasing Reynolds number, 
the j factor was found to rise to a maximum and then taper off. Figure 


21 indicates that -this peak in the j factor occurs between Reynolds 
numbers (based on Kb) of 0.5 and 1.5 
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There are three possible explanations for this anomalous behavior: 
1. Experimental error 
2. Longitudinal conduction effects 


3. Thermal boundary layer does not behave in the 
same fashion as the hydrodynamic boundary layer. 


If the cause were experimental ut then it would have to be a 
systematic کے‎ as all the cores behave in a similar manner. Trost 
[30] observed the same sort of behavior for Solar No. 1, Solar No. 6, 
Stainless Steel Plate-Fin reference matrix, and the entire series of 
brass cores (with varying L/r,, ratio). He could offer no explanation, 
and none can be offered here from a systematic error viewpoint. It 
might also be noted that no such behavior (i.e., a peaking of the j 
curve) was observed 70 the testing of the Solar No. 4 core. This 
would seem to rule out systematic error as a cause. 

| The inclusion of non-zero longitudinal conduction would tend to 
shift the low flow rate portion of the j curve up to where it would 
be more nearly parallel to the f curve, as demonstrated by Ball [2]. 
It has already been established that an improper assumption of À = 0 
will cause substantial error in the heat transfer results at low flow 
e. However, the small cross-section to surface area ratio and the 
sinuous path of the fibers make it difficult to conceive of such a 
mechanism at work in fibrous materials. A simple experiment to test 
this ہہت‎ is discussed in the section on suggestions for further 
investigation. 

The third explanation is the most plausible. When the thermal 
boundary layer no' longer behaves analagously to the hydrodynamic 
boundary DES, the Reynolds analogy is clearly invalid. This fre- 


quently happens in.external flows, where the turbulence induced by 
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the body contributes to the drag on the body, but only the skin 
former lon 1s ی‎ for the transfer of heat. Significantly, 
Colburn [6] noted that "the use of a friction line for predicting 
heat transfer would be very unsafe for flow across tubes." In par- 
ticular, when separation oct: it is well known that the j and f 
factors no longer behave in the same manner [18]. Inasmuch as effects 
similar to separation have been postulated as the cause of non-Darcy 
flow [3, 9, 12, 26, & 27] and inasmuch as virtually all the data runs 
were in the non-Darcy sea. the occurrence of separation or similar 
phenomenom would seem to be the most probable explanation for the 
failure of the f and j curves to behave in a parallel fashion. 

Figure 21 was an -attempt to see how well the j factor for the 
four samples was correlated with a Reynolds number based on Kb. There 
is really too much scatter in the data to draw any firm conclusions. 
The j factor, like the friction factor, is inversely proportional to 
the heat transfer area, A. By the same reasoning as before, if A 
were more precisely known, the curves might be more closely grouped. 
It is of interest to note that Trost [30] experienced e same inability 


to correlate j with N, (based on a well-defined channel hydraulic 


Re 


radius) for six different brass cores. 
Figure 21 also shows that the optimum Reynolds number for steel 
wool as a heat exchanger (neglecting friction considerations) is in 


the vicinity of N "I l, since the j factor is a maximum. Therefore 


R 


more heat is transferred per pound of air flow in that range of Reynolds 


numbers, hence the term optimum N, for heat exchange. 


Re 
A figure of merit in evaluating heat exchanger surfaces is the heat 


transfer power at a given friction power level (all evaluated at standard 
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conditions of 500°F and 1 atmosphere), as in Fig. 22. The higher a 
curve is on such a plot, the better the surface it represents is as 

a heat exchanger. Under the conditions of an assumed average equi- 
lateral cross-section for each fiber (from which the total heat 
transfer area is computed), and the assumption of operation at stan- 
dard ns. Fig. 20 indicates that Grade 3 steel wool is the best 
surface of those tested. It appears that the higher the grade, the 
better it is as a heat exchanger. All grades of steel wool tested 
were evidently better heat exchangers (on a per cee ee basis) than 


was Solar No. 4, which was also plotted. 
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VII. EXPERIMENTAL UNCERTAINTIES 


No attempt was made to evaluate errors arising from failure to 
satisfy the idealizations discussed in Section IV. 

The remaining sources of error are uncertainty in physical con- 
stants, inaccuracy in geometrical measurements, instrumentation error, 
and error in the permeability XA. The method of Kline and 
McClintock [14] was used to determine the uncertainty in the final 


results. 


A. UNCERTAINTY IN PHYSICAL CONSTANTS 

Values for the physical constants of plain carbon steel were 
obtained from [18], while those for air were from [10]. The uncer- 
tainties (as defined by [14] and different from percentage error) in 


these values appear to be [23]: 


es + 0.5 


من 


OMS‏ کے 
se ser:‏ 
o‏ . 
Nn? 2.0 |‏ 
ko: + 0.5%‏ 


= 
wn 
oe 


prt 


The thermal conductivity, Kos appears only in the longitudinal conduc- 
tion parameter, , which in this study is assumed equal to zero. Since 
the material of construction is specified only as plain carbon steel, 
there is a fairly wide latitude in the value IR and in particular, 


c.- The latter may be a major cause of inaccuracy [19]. 
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B. GEOMETRICAL MEASUREMENT INACCURACIES 

Contrary to previous studies, this area is a major source of error 
due to the deformable nature of fibrous media and the randomness of the 
internal geometry. The porosity of the sample seems to have the least 
uncertainty, since it was determined by two independent (though approx- 
imate) methods (see Section III) which agree to within one percent. There 
is variation in the length of the matrix (L) across its cross-section. 
This length also varies with flow rate because a certain amount of de- 
formation occurs as the vacuum downstream grows larger. The free flow 
area, E depends on the approach taken to the Reynolds number and the 
uncertainty in it is therefore a function of the accuracy of the approach. 
The crossed-rod matrices model yields a value of A, approximately 3$ less 
than that computed by the other two methods. The heat transfer area, A, 
has the largest uncertainty due to its extreme sensitivity to the assumed 
fiber cross-section. There is an uncertainty of 25% just in the selec- 
tion of a mean fiber width. The assumption of 30% uncertainty in A may 
even be optimistic. However, as noted by London [20] errors in such 
geometrical measurements do not affect the numerical re of the N. a 
results. The Colburn j factor and convective heat transfer coefficient, 
h, are affected, but the same results may be obtained if future work on 
geometrically similar surfaces is based on the same geometrical assump- 
tions. 

The uncertainties assumed for the various geometrical measurements 


dress 
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من 


A : + 5.0 
A : +6 


W : negligible 


(W. is the weight of the matrix). 


C. INSTRUMENTATION ERROR 

Instrumentation errors were primarily due to the pressure measure- 
ments. The uncertainties suggested by Piersall [23] for the various 
pressure measurements were not justified based on experience in this 
study of fibrous media. The following uncertainties are suggested in 


ivewmon Prersall's values: 


Po: 2.0% 
OP: 220% 
QP: + 2.5% 

P 00۷۳ ۰ھ‎ 


Ba negligible. 


The orifice temperature was measured with a sheathed copper- 
constantan thermocouple and read by a Leeds and Northrup potentiometer. 
Based on the manufacturer's calibration of the thermocouple, the esti- 
mated error is + one-half of the smallest potentiometer gradation, or 
5۳5 ۰ ۵ 02۰ ۱۱۱۱ 6۵6 (approximately + 0.1°F). The transient tem- 
perature response (t,-t,) was measured directly from the recorder out- 
put in millimeters. The assumed uncertainty in the maximum slope thus 
measured is + 3.0%, 

Since the N „y determined is a function solely of the maximum slope 
(for A= 0), it is appropriate to evaluate the uncertainty at this point. 


The bulk of the N. Values were in the range 22?N, 75. From Fig. 3 
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for A = 0 and an uncertainty in the maximum slope of 3%, the uncertainty 

in the N, is between 3.39 (N, = 22) and 3.75% (N. = 5). A few values 
Lu tu Tu 

of 8 were in the vicinity of 2.0, where the uncertainty in the Na is 


indeterminate. These results must therefore be regarded with skepticism. 


D. PERMEABILITY ANALYSIS ERROR 

Although primary data was used in the determination of the permea- 
iue. K, tortuosity meter. b, and cubic coefficient, e, of each 
TR the uncertainty limits cannot be found by the method of Kline 
and McClintock, since these quantities are not the result of single or 
multiple sample experiments, but rather a least squares curve fitting. 
The uncertainties in K, b, and c are a function of the spread (i.e., 
uniform distribution 7 the mass flow range) and number of data points, 
as well as the accuracy of the individual points. For these reasons 


the relative uncertainties in K, b, and c at 20:1 odds must be estimated: 


K: RN 
b: * 10$ 
C + 26% 


Ee RESULTS 


The basic formula for the uncertainty w, in result RIV, > V 


R ae 


"DE 
at specified odds is: 
! 
38 (AR, mp 
W, = F + CN ee ( We | 66 
R 3 ۷ ۱ aV, IVA (66) 
where Wi» WS, i... W. are the uncertainty intervals at the same specified 


ódds for the n variables of which R is a function. Note that Wp has the 


units of R. If the n variables appear in some form in the formula for 


7 





R as factors, rather than addends, simplification in (66) results 
from dividing through by R. 

Consider, as an example, the evaluation of the uncertainty in the 
mass flow ee, 7. (See Appendix B for the origin of the following 
equation): 


2 oU 
m = 0 (67) 


Let P. = Pse - P ° With (66): 





2 2 
TOM =(const)-| ES A + A Fo & ge 


۶ 
ANS (68) 
+AR E (-4 کے‎ ۱ 





Divide through (68) by (67): 


Win | 1 [eY ES E 
m +5 4 ۰ ۳ ZN 2 


Wwe 








The following relative uncertainties (i.e., 3 ) have already 
been postulated (based on 20:1 odds): 

= - e = 0.02 
Py m 
W, 

(=) = 0.02 
AP, 

(=) OIDO (negligible) 
Em VM 

0 O 


Therefore the relative uncertainty in the mass flow rate is: 


e - +14), 


In a similar manner, the relative error in the mass flux G is found to 


be 7 
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There is no reason to compute the uncertainty in the Reynolds number 
based on either the wa ٦ porosity or crossed-rod matrices 
models, as these are completely artificial concepts. Conversely, the 
Reynolds number based on the product Kb does have physical significance 
۳ there is a correct value). The uncertainty for this Reynolds 


number is: 


The uncertainties in the other results are as follows: 


11.4% 


Q 
{+ 


30.4% 


| + 


J Q am o. 

It should be noted once again that these uncertainties are based on 
20:1 odds and are not maximum percentages of error. No attempt has been 
made to include the effect of such theoretical assumptions as zero 
longitudinal conduction, step change in temperature, etc., in the 


uncertainty analysis. 


99 





VIII. CONCLUSIONS 


Based on the experimental results of this study, the following 


conclusions may be made: 


de 


The pressure gradient in fibrous media may be represented by 
a cubic equation with the bulk velocity as the independent 
variable: 


d = U, bu, "لا‎ 
EIKE de 


K is the permeability of the medium, b is a tortuosity 
parameter, and c is a small, negative, dimensionless number 
whose physical significance is unknown. The above equation 
assumes negligible compressibility effects, as found in 
this study. 


The product of the permeability and the tortuosity parameter 
is superior to the artificial hydraulic diameter determined 
from either the uni-dimensional porosity or crossed-rod 
matrices models for use as the characteristic length in 

the Reynolds number for flow in fibrous media. The product 
of the permeability and tortuosity parameter is a function 
only of the medium, but is not a conceivable dimension as 

is a hydraulic diameter. Moreover, this quasi-length parameter 
must be determined experimentally, and is very sensitive to 
the subleties in the structure of the fibrous medium. In 
light of the above conclusion that there is a cubic term in 
the equation for the pressure gradient, it is possible that 
the precise characteristic length may involve the dimension- 
less number c as well. 


A fiction factor correlation, as a function of N 2 based on 
Kb, which is valid for steel wools at Nee less than unity is: 


E - 207 3 Z, 


Por Noe in excess of unity, this relation over-estimates C 


E" 
The Colburn j factor shows a peak value in the vicinity of 

Reynolds numbers of one for all grades of steel wool. This 
is the optimum operating point in terms of heat transferred 


per pound of air flow. 


The thermal and hydrodynamic boundary layers in. fibrous media 
do not behave in analagous manners. Separation, which causes 
deviation from the Darcy flow regime is the probable explana- 
tion for this dissimilar behavior. 
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Grade 3 steel wool (by far the largest fiber size) is the 
best overall heat exchanger, based on evaluation at standard 
conditions and an assumed equilateral cross-section. In 


Memes the acces the fiber size, the better the perfor- 


mance as a heat exchanger. 
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IX. RECOMMENDATIONS FOR FURTHER STUDY 


This first-cut analysis of heat transfer and flow friction charac- 


teristics of fibrous media has raised questions which merit further 


investigation. 


1. 


' Since the surface area, A, was the largest source of 


uncertainty in the experimental results, the use of gas 
diffusion and other more sophisticated techniques for 
determining specific surface [27] should be pursued. 


The possibility of longitudinal conduction should be 
investigated. Since Grades 1, 2, 3 may be either stain- 
less steel or plain carbon steel, and since these two 
metals differ greatly only in thermal conductivities, 

a simple experiment to determine if any longitudinal 
conduction exists is easily performed. 


Inasmuch as values of N, less than 3.0 were experienced, 
the use of the cyclic testing technique for future studies 
is recommended. 


In addition to sampling various grades of steel wool, it 
is suggested that the same grade of steel wool be tested 
at various compactnesses. Thus porosity, as well as fiber 
size becomes an independent variable. 


Some further thought must be given to packaging a steel 
wool sample before testing, so that better control over 
homogeneity, dimensions, and boundary conditions may be 
had. 


Prior to any further use of the single-blow transient testing 


apparatus, some improvements in the instrumentation should be made. 


ir 


Radiation shielding should be applied to all thermocouple 


‚and heater leads to eliminate noise problems. 


Water manometers should be replaced with pressure transducers. 
In this way a graphical output is obtained, eliminating the 
need to make a value judgement under the press of time (i.e., 
during the course of a test). 


Additional pressure transducers should be installed down- 
stream of both the matrix and the orifice. These would 
Serve as a check, since the inlet pressure plus the 
pressure drop should equal the outlet pressure. 
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If all of the above instrumentation improvements are made, 
six pressures would be measured. Since an eight channel 
recorder is available, the orifice temperature could be 
placed on the seventh channel. Thus a complete graphical 
display for each data run would be available. 
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Figure lA.  Zxperimental Setup 





Figure 43. Detail of Inlet Cone, ‘Jeavers, Test Section 
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Fanning Friction Factor as a Function of Reynolds Number 


Based on the Crossed-Rod Matrices Model. 


Figure 10. 
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Fanning Friction Factor as a Function of Reynolds Number 


Based on the Product of the Medium Permeability and 


Tortuosity Parameter (Kb). 


Figure 1l. 
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Grade 00 Steel Wool 
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Figure 13. Geometrical and Physical Properties of 


Grade OO Steel Wool 
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Figure 1h. Heat Transfer and Flow Friction Characteristics of 
Grade OO Steel “lool 
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reometrical and Physical Properties of 











Figure 16. Heat Transfer and Flow Friction Characteristics of 
Grade O Steel Wool 
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(Divisions in mm. 


Grade 1 Steel ‘lool 
FRONTAL AREA FH 2 uÁ 01056 
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Figure 17. Geometrical and Physical Properties of 
Grade 1 Steel “ool 
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Figure 18, Heat Transfer and Flow Friction Characteristics of 
Grade 1 Steel Wool 
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Figure 19: Geometrical and Physical Properties of 


Grade 3 Steel wool 
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Figure 20, Heat Transfer and Flow Friction Characteristics of 
Grade 3 Steel Wool 
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Figure 21. Colburn j-Factor Versus Reynolds Number 
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APPENDIX A: DESCRIPTION OF EQUIPMENT 


A. AIR SUPPLY AND FLOW MEASUREMENT 

A 30 HP, multi-stage Spencer Turbo-Compressor rated at 550 CIM 
(for intake at atmospheric pressure) provided an induced draft through 
the test section. 

The air flow rate was measured with an ASME standard orifice section. 
The concentric orifices had throat diameters of 0.308, 0.462, 0.775, 
10232, 1.540, AS 7I and 2.310 inches. These values correspond to 
throat diameter to pipe diameter ratios B) 657۶0۰۲0۱۰۸۰۰۷۰۹ 77ء‎ 
0.40, 0.50, 0.64, and 0.75 respectively. Pressure taps were located 
one pipe diameter upstream and one-half pipe diameter downstream of 


the orifice plate, as specified by [1]. 


B: HEATER SYSTEM 

The heater system used was a that used by Bruno [5], 
although it was completely rewired. Fourteen separate bakelite frames 
were each wound with two parallel-connected heater elements. The heater 
elements were 0.0031 inch nichrome wire, (Fig. 23) which was selected 
largely due to the comparatively low coefficient of thermal expansion 
of nichrome. Each heater element has a nominal resistance of 800 ohms. 
An individual switch controlled both heaters per bakelite frame (i.e., 
14 such switches) and these switches were in turn controlled by a master 
switch which was used to cut off the power to the heaters to simulate 
the step change in air inlet temperature . Thus the number of heaters 
(as well as the voltage supplied) could be adjusted to provide the 20°F 
temperature rise for any flow rate up to approximately 1,000 pounds per 


hour. 
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C. HEATER POWER SUPPLY 

A General Electric Speed Variator motor-generator, operating on 220 
volts AC at 35 amps and producing 250 volts DC at up to 28 amps, provided 
the electrical power to the heaters. The actual voltage supplied was 
regulated with the rheostat on the field winding of the MG an and at 
the test section with a 100 ohm Ohmite wound resistor (for currents less 
than 4 amperes). 

The current and voltage at the source were monitored with a Roller- 
Smith 150 a. ammeter and a Simpson multimeter, respectively. At the 
test section the current and voltage were Ber cone with a Westinghouse 


10 a. ammeter and 300 v. voltmeter. 


D. MATRIX HOLDER AND TEST SECTION 

Both the holder and test section (Fig. 24) were fabricated from 
plexiglass. The matrix mde is a drawer which slides into the test 
section. By removing the rear, ‘upstream, and downstream faces of the 
holder, the test section of steel wool may be inserted with minimum 
packing effort, thereby causing negligible inhomogeneity in the porosity. 
The flow channel is 3.1875 inches square and samples up to approximately 
3 inches long can be handled. A network of seven .003 inch steel wires 
was strung across the downstream face of the holder to restrain the 
steel wool from being sucked down the duct. These wires were not thought 
to interfere with the friction or heat transfer performance of the sample, 
and experience indicated that they are probably unnecessary anyhow. A 
Similar network was placed in e upstream face to prevent interference 
between strands of steel wool and the upstream thermocouples. This net- 
work is also thought to be unnecessary. Styrofoam insulation lined the 


interior of the matrix holder and also served to hold the -steciewe olean 


place. 
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The test section, into which the holder slides, contains the matrix 
upstream and downstream thermocouples and کت‎ taps. In previous 
installations the matrix downstream temperature was measured by thermo- 
couples mounted in the as rather than the test section. However, 
the possibility of shorts or breakage from pertruding strands of steel 
wool dictated the movement of the thermocouples downstream by approxi- 
mately 0.5 ine No significant error in the maximum slope technique 
is thought to have been introduced by this movement of the most impor- 


tant sensor. 


E. INLET CONE AND FLOW STRAIGHTENER 

This section was designed and tested by Piersall [23], and was 
found to provide a uniform velocity profile for the air entering the 
matrix. 

Additionally an alteration was made in the radiation shielding of 
the reference thermocouples, located at exit from the inlet cone, which 
may have improved the uniformity of the velocity profile. Originally 
all the reference thermocouples were contained in one 3/8 inch aluminum 
tube mounted across the center of the channel. In this reconstruction 
these thermocouples were placed in five 1/8 inch nickel-plated tubes, 
evenly spaced across the channel. This is thought to delocalize the 
intensity of the disturbance caused by the presence of the thermocouples 


in the stream. 


F. PRESSURE MEASURING SYSTEM 
One pressure tap was located in both the upstream and downstream 
E o ol Ile test ction (i-e. on either side of the matrix). As 


previously mentioned, two taps were also located at the orifice section. 
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Each pressure tap was connected by Poly-flo flexible tubing to a high 
pressure or low pressure manifold and 00 by a Hoke quick-closing 
valve. The high and low pressure manifolds were connected to the high 
and low pressure sides of the following instruments: 


2 Merriam Instrument Company, 50 in. manometers 
1 Ellison Draft Gage Company, 0-2 in. inclined manometer 


1 Ellison Draft Gage Company, 0-4 in. inclined manometer. 
These manometers were also protected with quick-closing valves. Thus, 
with the proper valve Er each of the four pressures measured per 
data run was simultaneously measured on four manometers until the mass 
flow rate and orifice plate combination required that the more limited 
manometers be isolated. 


Atmospheric pressure was measured by a PRINCO mercury barometer. 


G. TEMPERATURE MEASURING SYSTEM 
Temperature sensors were located at the exit from the inlet cone 


) و(‎ entrance to the test section i c exit from the test section 


(t4), and 15.5 inches upstream from the orifice section (t ). 


was measured by two different sets of five 30 gage 


Temperature t 


iron-constantan thermocouples connected in series. Each thermocouple 
was wrapped in teflon tape to prevent shorting. As previously mentioned, 
the reconstructed apparatus has different housing for these thermo- 
couples than that designed by Traister [29]. The ten thermocouples 
were housed two to each of five nickel-plated brass tubes. The 
two thermocouples in each tube were from the two different sets of t 
sensors. The tubes were cut-away to the front while closed to the 
downstream side, where they faced the heaters. This system is thought 


to provide superior sensitivity and radiation shielding to the old systen, 


where all ten thermocouples were housed in one aluminum tube. 
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Temperature t, was measured by a set of five .001 inch diameter 


2 
iron-constantan thermocouples in series. The output of this thermopile 
was referenced against one of the t; thermopiles so that the net output 
to the recorder was (t,7t,). 


Temperature t, was measured exactly as t so the output to the 


3 28 


recorder was (t, -t, ). 

The two outputs (t,-t,) and (t,-t,) were amplified by Hewlett- 
Packard 8803A DC amplifiers and recorded on a Hewlett-Packard 7702B 
two-channel recorder. 

The orifice temperature, er was measured by a 30 gage cooper- 
constantan thermocouple sheathed in stainless steel. It was referenced 


to an ice junction and the output read on a Leeds and Northrup Milli- 


volt Potentiometer. 


H. NOISE CONTROL SYSTEM 

A substantial amount of noise in the (t,-t,) and (t,-t,) signals 
was frequently encountered. To filter this out, Mallory 2.0 microfarad 
capacitors were placed between the negative side of the signal and 
ground. The time constant of the circuit with the addition of these 
capacitors was less than the minimum time constant of the heater-thermo- 
couple system reported by Bruno [5]. 

Additionally a Sprague 1000 micro-farad capacitor was placed in the 
power supply circuit to filter any noise from the motor-generator. This 
had no effect on the transient response of the apparatus, since the 


master switch isolated the power supply from the heaters. 


99 





APPENDIX B: DATA REDUCTION RELATIONSHIPS 


A. GEOMETRICAL RELATIONSHIPS 

Once an average cross-section for a fiber is assumed, the heat 
transfer area of the medium is determined. In this study the average 
cross-section was assumed to be an tral triangle with side 
length w, where w is the median mean width of the fiber specified 
for the particular grade steel wool by [7]. Therefore: 


Cross-sectional area = A.. - w*(¥3/4 ) 
fiber 


Perimeter of cross-section = 3w 


Aé 
Hydraulic diameter of fiber = ج‫ جآ‎ = Ww 3 (B-1) 


Volume of steel in sample - V er on) (Total length of fibers) 


steel” fib 


Heat transfer area = A = (Perimeter)(Total length of fibers) 


(Sw) ات‎ /A ç. m E NS Va In 


But the volume of steel is equal to the weight of the sample divided by 
the density of steel. te. 
Vsieel ü Ws /Ps 
A ٤ Ws ا یع‎ 
Porosity = ff = Ne qus Ws Iles AsL) (B-3) 


Specific surface - B= AÑ. = ۸۷۷3 W. Kg Aj) (B-4) 


Free flow area = A ^ (except for crossed rod matrices) (B-5) 
f € fr 


(B-2) 


Equations (B-1, B-2, B-3, & B-5) are substituted into equations (1) through 
(6) as appropriate to determine the geometrical characteristics for use 


in the various models. 
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B. MASS FLOW RATE 
The mass flow rate was calculated from the ASME Power Test Code [1] 


with Murdock's [21] interpolation scheme. 
š = s he Y 
و‎ 
۳۳۱ 3 ol = E EY o (B-6) 


where 


C = orifice coefficient of discharge 
6 = ratio of orifice diameter to pipe diameter 
d = orifice diameter in inches 
F_ = orifice plate thermal expansion factor 
= fluid thermal expansion factor 
= specific weight of fluid flowing = P/RT if a perfect gas 
P = absolute static pressure at orifice (3b£/ft^) 
R = specific gas constant for air = 53.35 ft-1bf/1bm-"F 
T = absolute temperature at ee (AR) 


AP = pressure drop across orifice in inches of water 


Substitution of the identities for K and in (B-6) yields: 


T‏ سا rd‏ سی 


From"Fig. HOA ot [l1] one finds: 
A S 
pae der O 350 Vic 


where 
ks c,/c, > ۱۹4 (for air) 


x - AB/P (any consistent units). 


JO] 





Pico from [i]. Fig. 38: 


۳  < 1 
3 
= (P , - P_/13.6)(70.75)1bf/ft” 
atm O 
TES local atmospheric pressure in inches of Hg 
P = static pressure upstream of orifice plate in inches H,0 


2 


When all these identities are substituted into (B-7) with the units 


adjusted to give m in (lbm/hr), the equation becomes: 


۱ ۱ 
382,810 [42-041 + 0358] AP |aR(Aim- Balız.c)]” 


Tarp | 990 (Rin Pao) | to + 459.7 


The orifice coefficient of discharze is found from Murdock!s tables 


and interpolation scheme. 


C. MAXIMUM SLOPE 
The following sketch of the general time-temperature history down- 


stream of the core will illustrate the relationships involved. 


/ 
4 


c 
۱ 
Ct 
ST 
e 
| 
3 
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The maximum slope of the general time-temperature history is a function 


of N and A. 
u 


u) 
AN) 


MAY 





= cp NR) 


The non-dimensional variables are: 


T = generalized time variable = LAG /(Ws¢s\ 
N ° hA/me, 


TIN z [eoe YO 


Therefore 


: ۰ _ نم 


Furthermore 
t; - t. E t.-t, En 1 
to- t. ID RUE 


and its derivative is 


UN t; k 4L dá 
k xus 


t 


Combination of (B-9) and (B-10) yields: 





C 
ag zm ME d(t,-t.) 


uam) me 48 
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(B-9) 


(B-10) 


(BS) 





In terms of the recorder output the various quantities in (B-11) are: 


Elan] = Y/X 


x/(Chart speed) = 0ے ت7‎ 


d(t,—t,) Y 


II 
3 


(t,-t,) 


Equation (B-11) becomes in these terms: 


EZ pe 
E ela} (B-12) 
MAX 


This value of the maximum slope and the known value of A(in this study, 
A= 0) are used to enter Table II or Fig. 1 to determine the correspond- 


ing value of Nea 


D. COLBURN J FACTOR 


The Colburn j factor is defined as 


2/5 h 2/3 
= N N = — 
J st! ‘Pe Gc, No, (B-13) 


In this case the mass flux, G, is based on Du rather than A ep’ although 
at 98% porosity there is not much difference. Substitution of G = m/A 


and multiplication of (B-13) by A/A yields: 


ETS 
J~ me, A Ne. 





But N =" (Arc ) and therefore: : 
tu p 


Š No, (=) NE‏ ل 
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E. HEAT TRANSFER POWER AND PRICTION POWER 
The heat transfer power per unit area per degree temperature 


difference is: 





hz Eo (Na | (B-14) 


Pe 


(where 1 is the characteristic length used in PF 


At standard conditions of 500°F and one atmosphere [13] 


0.2477 BTU/1bm-°F 


0.0678 lbm/hr-ft 


9 OnT کر‎ 


<p 
مر‎ 
j 


N 


Pr 0.671 


Substitution of these values into (B-14) leads to: 


BTU 


-0.02195(4\N, 7 لا‎ 
hesg = O.OZI9S ( N | zZ Dee (B-15) 


The flow friction power per unit of area is [13] 
5 3 3 
5. (T) 2 N 
E- faz ] >) Re. 
Evaluation of (B-16) at standard conditions yields: 


T A "T ۳ 2 
2-07 "03 [Q F |o? "n (B-17) 


(B-16) 





P. OTHER RELATIONSHIPS 
The relationships pertaining to the Reynolds number and the two 
friction factors have been completely derived in the section on 


theoretical background. 
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